Quantum discord in finite XY chains 
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We examine the quantum discord between two spins in the exact ground state of finite spin 1/2 
arrays with anisotropic XY couplings in a transverse field B. ft is shown that in the vicinity of the 
factorizing field B s , the discord approaches a common finite non-negligible limit which is independent 
of the pair separation and the coupling range. An analytic expression of this limit is provided. The 
discord of a mixture of aligned pairs in two different directions, crucial for the previous results, 
is analyzed in detail, including the evaluation of coherence effects, relevant in small samples and 
responsible for a parity splitting at B s . Exact results for finite chains with first neighbor and full 
range couplings and their interpretation in terms of such mixtures are provided. 

PACS numbers: 03.65. Ud,03.67.Mn,75. 10. Jm 



I. INTRODUCTION 

The great interest on quantum entanglement in recent 
years was triggered by the key role it played in certain 
quantum processingtasks like quantum teleportation and 
superdcnsc coding [13]. It was also shown to be essen- 
tial for achieving an exponential speedup over classical 
computation in pure-state based quantum computation 
However, it was recently recognized that in the case 
of mixed-state based quantum computation, like the "de- 
terministic quantum computation with one qubit" model 
introduced by Knill and Laflammc Q, an exponential 
speed up can take place without a substantial presence 
of entanglement @. This has turned the attention to 
alternative measures of quantum correlations in mixed 
states, like the quantum discord, introduced by Ollivier 
and Zurek Q, which can detect those quantum corre- 
lations present in certain separable mixed states 0] and 
hence not captured by the entanglement of formation Q , 
but still useful and crucial for certain quantum tasks. It 
was in fact recently shown by Datta, Shaji and Caves [Io| 
that the circuit of Q does exhibit a finite non-negligible 
value of the quantum discord between the control qubit 
and the remaining mixed qubits. Since then, interest on 
the quantum discord and other alternative measures has 
grown considerably pH4l5| and several studies of their 
behavior in spin p airs immersed in a spin chain have been 
performed [16l4l9l |. 

The aim of this work is to analyze the quantum discord 
of spin pairs in the exact ground state of finite spin 1/2 
arrays with anisotropic XY or XY Z type couplings in a 
uniform transverse magnetic field B. The exact ground 
state of a finite chain will have a definite spin parity and 
this property will be seen to deeply affect the discord for 
fields lower than the critical field B c , where we will show 
that the main results can be interpreted in terms of the 
discord of mixtures of aligned pairs. Our interest in these 
systems is motivated in particular by the remarkable fac- 
torization phenomenon they can exhibit at a particular 
finite value B s < B c of the magnetic field [20 - 1271 ]: At 
such field, they have an exactly separable (i.e., factor- 
ized) ground state. This feature was first discovered in 



[20] in one-dimensional (1-d) XY chains with first neigh- 
bor couplings, and later shown to occur also in more gen- 
eral systems, like 2-d arrays [22|], cyclic 24 1 and general 
[251 . 12a ] arrays with arbitrary range couplings with a com- 
mon anisotropy and also systems in non-uniform fields 
[2^ . For transverse fields, these separable ground states 
actually break parity symmetry and are hence degener- 
ate, coinciding B s in a finite system with the last cross- 
ing of the two lowest opposite parity levels (24|. A most 
remarkable related feature is that in the immediate vicin- 
ity of B s , pairwise entanglement, though weak, reaches 
full range [H, [24[ , regardless of the coupling range [24| , 
changing at B s from antiparallcl (B < B s ) to parallel 
(B > B s ) type (23[, indicating an "entanglement transi- 
tion" . This suggests the possibility of a non-zero discord 
between distant pairs at least in the vicinity of B s , with 
universal features (independence of separation and cou- 
pling range). Here we will show that this is indeed the 
case, and derive also its analytic limits at this field. 

Moreover, distributed pairwise entanglement is neces- 
sarily weak due to the well-known monogamy property 
[28l.[29j associated with entanglement sharing: If in a sys- 
tem of n spins or qubits one spin is strongly entangled 
with a second spin, it cannot be strongly entangled with 
any of the remaining spins. This fundamental feature 
follows from the bound [H, H| Y^Cfj < C? < 1 satis- 
fied by the pairwise concurrences [30] measuring the 
entanglement between spins i and j, where Cj represents 
the concurrence between i and the rest of the chain. If 
all pairs are equally entangled, the maximum value that 
can be reached by CV,- is in fact just 2/n [HI, [Hj]. In 
contrast, the quantum discord is not affected by such 
bound and can simultaneously reach non-negligible finite 
values between any two spins, as will be seen to occur in 
the vicinity of B s , leading to a quite different behavior 
with the applied field and separation in the whole region 
|£>| < B c . The properties of the ground state discord 
in the vicinity of the separability field in anisotropic XY 
chains were not discussed in previous references. We will 
also analyze here finite size effects, which lead to a finite 
step of the discord at the factorizing field B s and other 
parity transitions, visible in small chains. 
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Sec. |TT] discusses the quantum discord and its evalua- 
tion in typical reduced states of a spin pair in such chains, 
describing in detail the case of a mixture of two aligned 
states, which represents the exact reduced state in the 
vicinity of B s . Coherence effects in these mixtures, rel- 
evant for small chains, are also examined. Sec. IIIII dis- 
cusses the behavior of the discord of spin pairs with the 
applied field and separation in finite cyclic XY chains 
with first neighbor as well as full range couplings, in- 
cluding its interpretation and the differences with the 
pair entanglement. The appendix discusses the details of 
the exact definite parity solution of the finite cyclic XY 
chain. Conclusions are finally drawn in IIV1 

II. FORMALISM 

A. Quantum Discord 

The quantum discord (D) is a measure of quantum 
correlations based on the difference between two distinct 
quantum generalizations of the classical mutual infor- 
mation, or equivalently, the classical conditional entropy 
@, EH • Given a bipartite system A + B in a mixed state 
Pab, and denoting with p A = Tr B p AB , Pb = T^aPab 
the reduced density operators of each subsystem, D can 
be expressed as Q 

D = I{A:B)-M&xI {p b } (A:B) (la) 

= Mm[S(p' AB ) - S(p' B )} - [S(pab) - S(p B )} ,(lb) 
{Pf} 

where I(A : B) = S(p A ) + S{p B ) - S(pab) > is the 
quantum mutual information, with S(p) = — Trplog 2 p 
the von Neumann entropy, and I^ p b^(A : B) = S(pa) + 

S(Pb) — S{p' AB ) the mutual information after a local von 
Neumann measurement in system B defined by a set of 
orthogonal projectors Pf = Ia® |js)0's|- Here 



represents the joint state after such measurement if the 
result is unknown, with p' B = Tta p'a B an d Pa = 
Tt b p' AB — pa- Minimization in (|lb[) is over all sets of 
local projectors. 

The last bracket in (jlbj) is the direct quantum 
extension of the classical conditional entropy 0, HH 
(S(A\B) = S(A,B) - S(B) = EjPjSiA/j), with 
S(A,B) = Y^p.Ao^p,,. S(B) = -J2 jP flog 2 pf 

and S(A/j) = -J2j(Pij/pf) lo &2(Pij/pf) for a system 
with joint probability distribution pij and marginal dis- 
tribution p® = 'Y^ li Pij) 1 which, in contrast with the clas- 
sical case, is not necessarily positive. On the other hand, 
the first bracket is the conditional entropy S^ p b^(A\B) = 
PjS(pj) after a local measurement in system B, where 
Pj = Tr P^ pab is the probability of outcome j and 



Pj = P^ pAsPf ' IPj the state after such outcome. It rep- 
resents the average lack of information about A after such 
measurement and is a non-negative quantity. The second 
term in (|la[) is considered a measure of the classical cor- 
relations between A and B in [34| , the discord measuring 
then the quantum part of these correlations. 

In the case of a pure state {p\ B = pab), S(pab) = 
while the first bracket in (|lb|) vanishes for any choice of 
local projectors, and the discord reduces then to the en- 
tanglement entropy [35| D = E = S(pa) = S(ps)- In the 
case of mixed states, however, it does not coincide in gen- 
eral with the entanglement of formation Q (the convex 
roof extension of the entanglement entropy). While the 
latter vanishes for any separable state, i.e., for a convex 
superposition of product states pab — J2 a 1aPA ® P%, 
Qa > HI, the discord can be non-zero for these states. 
The discord vanishes in the case of densities diagonal in 
an orthogonal product basis (bAB = {Ka)|j_b)}) or in 
general a conditional product basis (bAB = 
with the orthogonal set {I^a)} depending on |js)), but 
will not vanish in general for a mixture of non-commuting 
product states, which is still a separable state. The dif- 
ference |T]) can be shown to be non-negative 0] due to 
the subtle concavity property of the conditional von Neu- 
mann entropy S(pab) ~ S(pb) [33^ 

B. Quantum discord and entanglement of spin 
pairs in definite parity states 

Let us now describe the basic elements to evaluate the 
discord and entanglement of a spin pair in a typical eigen- 
state of a finite chain of n spins, where the rest of the 
spins will play the role of an environment. We will con- 
sider spin 1/2 chains with XYZ couplings of arbitrary 
range in a uniform transverse magnetic field B along the 
z axis, such that the chain Hamiltonian has the form 

H = Bj2s iz ~hYl E J ^V^, (3) 

where Si M denotes the spin components at site i (in units 
of h). H commutes with the spin parity operator 

P z = exp[i7r J>« + V2)] = IJ(- <7 «) > ( 4 ) 

i i 

where o~i Z = 2si Z , which changes Si^ to — Si^ V i and 
p, = x, y. Non-degenerate eigenstates of H will then 
have a definite spin parity P z = ±1, a symmetry which 
will play a fundamental role in our discussion. 

The reduced density matrix of an arbitrary pair i, j in 
such eigenstate, 

Pij =Tr„_y|* v )(*„|, (5) 

will then contain no elements connecting states of op- 
posite parity, commuting therefore with the pair parity 
Pj* = a iz a jz . In the standard basis {|00), |01), |10), |11)} 
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(\ki) = \h)\i,), with aulki) = ^ ki \h), h = 0,1), Pij 

will then be of the form 



Pij 



a a 

c /3 

/3 c' 

a 6 



(6) 



where, setting (O) = trp.y O and s^± = s x ± is^, 
(?) = i±i(s. 

(^) = i±5( s 

{%) = • 



jz) + {SizS jz ) , , (7) 



^IZ°JZ/ ) 



(8) 
(9) 



with a + 6 + c + c' = 1. We will here consider translation- 
ally invariant systems such that (si Z ) = (sj Z ) and hence 



1(1 



b). Moreover, a and /? will be real 



since H can be represented by a real matrix in the stan- 
dard product basis of the full space. Non-negativity of 
Pij implies |a| < \fab, < c, with a,b,c non-negative. 

The internal entanglement of the pair can be measured 
through the entanglement of formation, which for the 
case of two qubits can be explicitly calculated as [13] 



E = -J2l^Og 2 q Ul q ± = \{l±.Jl-C 2 ), (10) 



where C is the concurrence [30j . It is given here by 



C = 2Maxf|a| 



ab,0] 



(11) 



The entanglement of the pair is of parallel (antiparallel) 
type if the first (second) entry is positive [23, [24| . Just 
one of these entries can be positive for a non- negative pij . 

On the other hand, in order to evaluate the discord of 
the pair, we need first the eigenvalues of pij, given for 

/ 1 — a — b 
\ 2 



d by A l 



± 



a+b 



± 



H] \ 2 If H 2 V V 2 

those of the single spin density matrix, 



^) 2 + |a| 2 ), and 



Pj = TT iPv 



a + c 
c + b 



(12) 



which are obviously Xj = |[1 ± (a — b)]. We also need 
to consider a measurement of the spin at site j along an 
arbitrary axis z' defined by the angles 7 and 0. The state 
of the pair after such measurement (Eq. ([2])) is 



p i i j=n,p ij pi,+pi.p ij pi, as) 

= U® \k'){k'\ for k = 0,1 with 

|0')=cos(i 7 )|0) + e ^sin(i 7 )|l) (14) 

|l')=cos(i 7 )|l)- e -^sin(i 7 )|0), (15) 



where Pi, 



\e™ k \k!) 



for Sj z i 



Sjz cos 7 



such that Sj Z '\k') 

Sj x sin 7 cos + Sj y sin 7 sin 0. For real a, (3, the eigen- 
values of pj; are 



A'; 



l+i/(a— b) cos , 7+A i x /[(2(a+b)-l)cos 7 +;y(a-b)] 2 +4|a+/3| 2 sin 2 7 

(16) 



where v = ±1, /1 = ±1 and |a + (3\ 2 = a 2 + (3 2 + 
2a/3cos2(j), corresponding v = 1 (—1) to the eigenval- 
ues of the first (second) term in (fT3"|). The eigenvalues of 
p'j = Tti p\j are then A^- = \ [1 ± (a — b) cos 7] . 

We have then all elements to evaluate the difference 

D( % 0) = S^j) - S{p'j) - [S(pij) - S(pj)} , (17) 

whose minimum (with respect to 9, 0) is the discord D 
(Eq. ([1])). For a/3 > 0, minimization with respect to 
yields cos 20=1 and just the minimization over 7 is 
finally required, which can be restricted to the interval 
[0,7r/2]. The minimum for the pair densities used in sec. 
Mil was obtained for 7 = ir/2 (z' ~ x), where = 

i+v\/( a >>) +i(a+W T_ ^, ecomes independent of v and hence 
degenerate. A general evaluation of the discord for states 
of the form (J6j> was recently provided fl4j . 



C. The case of a mixture of two aligned states 

A particular case of ([6]) of exceptional interest is that of 
a statistical mixture of two aligned states along arbitrary 
directions, not necessarily opposite, such that the local 
states involved are non-orthogonal. Choosing the z axis 
as the bisector of the angle between the two directions, 
we can write this state as 



-0)(-6-6\) 



a a 

a a 

a a 

a b 



Vb> - 4(1 ± COS* 
a = j sin 2 9 



(18) 



,(19) 



where \6) = exp^s^O) = cos \0\Q) + sin \9\l) and (H]) 
is again the standard basis representation. Eq. (jT5)) is the 
exact reduced state of any two spins in the immediate 
vicinity of the factorizing field (see sec. IIIip if coherence 
terms are neglected. It also provides the basic approxi- 
mate picture of the pair state in the region \B\ < B c . 

The state (fl"8)) is obviously separable Q, i.e., a convex 
combination of product states, and therefore its concur- 
rence and entanglement arc identically 0, as verified from 
Eq. (|11[) . However, its discord is positive for 9 £ (0, 7r/2), 
vanishing just for 9 = or 9 = tt/2: If 9 = 0, Eq. (fTHJ) 
becomes a pure product state (and hence D( 7 , 0) = V 
7,0) whereas if 9 = n/2, \9) and | — 9) are orthogonal 
and p{9) becomes diagonal in a product basis, with Eq. 
(fi"7|) vanishing then for 7 = n/2 and = (p[j = p^). 

In the general case, the eigenvalues of (fT9|) are \j = 
(^(1 ± cos 2 9), 0,0), with those of Pj{9) given by Xj = 
\{1 ± cos^), whereas the ensuing eigenvalues (| 16[) of p' irj 
become, for cos 20 = 1, 



./ \-\-v cos 9 cos 7+/i-^/[(l+i^ cos 9 cos j) 2 cos 2 #+sin 2 7 sin 4 9 
\j = 4 

(20) 

with those of p'j given by A^ = |(1 ± cos 9 cos 7). For 
9 e (0,7r/2], the minimum of (TTT|) is always attained at 
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1 + pJl - ±sin 2 29 
D= £{[2fc( Li 

f*=±l 

-K^^)-K^f^)]}, (2D 

where /i(a;) = —a; log 2 x. D is maximum at 9 — 9 m sa 
1.157r/4, where D sa 0.15. As seen in the bottom panel, 
while S(p' i j) — S(p'j) (the first bracket in (JH])) is an even 
function of 9 — 7r/4, S(pij) — S(pj) (the last bracket) is 
not, being maximum at 9 « 0.9l7r/4 and originating the 
deviation of 9 m from 7r/4. 

For 9 — > 0, D vanishes quadratically {D w -i^ 2 ) while 
for ^ f , D « - ±(§ - 0) 2 [log 2 (§ - 0) 2 + log 2 e - 2]. 



D. Effects of coherence term 



FIG. 1. (Color online) Top: The difference ([T7j) for the mix- 
ture (|18[) as a function of 7 and <j> at fixed 9 — n/4 (top) 
and as a function of 7 and 6 at = (bottom). D(p/,<$) is 
minimum at 7 = 7r/2 and = V 6 6 (0,7r/2]. 




FIG. 2. (Color online) Top: The quantum discord D = 
D(-k/2, 0) as a function of 9 for the state (|18[) . D is maximum 
at = # m ~ 1.15-71-/4. Bottom: The conditional entropies 
S(p'ij) — S(p'j) (upper curve) and S(ptj) — S(pj) (lower curve) 
as a function of 9, whose difference is the discord of the top 
panel. 



7 = 7r/2 (and cos 20 = 1), i.e., for a local measurement 
along the x axis, as seen in Fig.[TJ with D(^, 0) becoming 
quite fiat for small 9. Both 7 = 7r/2 and 7 = are 
stationary points of D, with a maximum. Fig. [5] depicts 
the minimum D = D(ir/2,0), which is the discord, as a 
function of 9, given explicitly by 



The reduced state of two spins at the factorizing field 
actually contains a small coherence term tx e(\99)(—9 — 
9\ + I - 9~9}(99\) (see sec.|nj), which leads to the state 



-6-6){-0- 



■(\99)(-9-9\ 



2(1 



( a a 

0/3/3 

0/3/3 

\a b 



(?) 



:(99\-9-9)) 

_ (l+e)(l±cosfl) : 
~~ 4(l+scos 2 9) 

(l±e)sin 2 8 
~ 4(l+e cos 2 9) 



(23) 



where |e| < 1. This term generates then a parity de- 
pendent correction to the previous results. The eigen- 
values of pij and pj are now Ay = (a + 6,2/3,0,0), 
Xj = (a + /3,b + (5), whereas those of p'^ and p'j can be 
obtained from Eq. (j!6p . The minimum of (|17[) is again 
obtained at 7 = 7r/2 (and <p — 0), (the surface being 
again similar to that of Fig. [IJ, leading to 



(1 + e) 2 + sin 4 



4(1 + ecos 2 



(1 + pcos 2 9)(l + pe) (1 +pcos9){l + pecos9) 
KH on I TZZ?TF\ ) _ M on , 2 m Qw 



2(1 + £ cos 2 9) 



2(1 + £ cos 2 



For £/0a nonzero entanglement of the pair also arises, 
with concurrence 



C = 



£ sm 



1 + e cos 2 1 



(25) 



which is parallel (antiparallcl) for e > (< 0). 

In the limit £ —¥ ±1, Eq. (|22|) becomes a pure state, 
namely, — > |^ r ±)(*±| with 



I* 



± 



^2(1 ± cos 2 



^/(1+cos 2 0)/2 
|01) + |10) 
V2 



5 




" ■ , , i > T -> i ,—---T', , , " 

-1 -0.5 0.5 1 



FIG. 3. (Color online) Effect of coherence term. The quantum 
discord (D) and entanglement of formation (E) of the state 
(122tl as a function of 8 for e = 0.2, and —0.2 (top) and as a 
function of s at 6 — n/4 (bottom). Entanglement vanishes at 
e — but becomes larger than the discord close to the pure 
state limit e = ±1, where D and E coincide. 

Therefore, D and E merge V 9 in this limit. Whereas 
l^-) is a Bell state independent of 9 (for 6^0), leading 
to D = E = C = 1, the entanglement of \^+) depends 
on 9 (with C = sin 2 9/(1 + cos 2 9) < 1), increasing with 
increasing 9 £ [0,n/2]. 

The response of D and E to the coherence term is 
shown in Fig. [3] For sufficiently small e, the correc- 
tion to D is linear in e for 9 not close to ir/2, with the 
discord increasing (decreasing) for e < (> 0), while 
at 9 = 7r/2 the correction is quadratic and positive (at 
9 = tt/2, D = 1- £ M=±1 /i(^) « i£ 2 log 2 e). En- 
tanglement, on the other hand, becomes finite as soon 
as |e| increases, becoming even larger than the discord 
for 9 close to tt/2 (where C = |e| and E cx — e 2 logs 2 for 
small e). As seen in the bottom panel, at an intermediate 
9 entanglement remains smaller than the discord just in 
an interval around e = 0, becoming slightly larger before 
reaching the pure limit e = ±1, where D and E coincide. 

III. QUANTUM DISCORD IN XY CHAINS 

We have now all the elements to evaluate and under- 
stand the behavior of the discord of an arbitrary spin pair 
in the ground state of a finite chain. We first consider a 
finite 1-d XY cyclic chain of n spins with first neighbor 
couplings, where J l J = and 

J% =6j,i±iJn, fJ> = x,V, (26) 

with n + 1 = 1. The exact solution for finite n can be 
obtained through the Jordan- Wigner fermionization [3?jj 
(see Appendix). The exact ground state will have a defi- 
nite (field dependent) spin parity and the reduced density 
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n=100 
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D 
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n=100 
L=2 

D 

E 

r^f^ro — i — 


— i h— i 

*^\^ \ 


n=100 
L=3 

D 

E 




r-f ' — 1 — 1 — 

n=100 



1 2 

B/J x 

FIG. 4. (Color online) Quantum discord (D) and entangle- 
ment (E) of spin pairs with separation L in the exact ground 
state of a cyclic chain of n — 100 spins with first-neighbor XY 
couplings and J y /J x — 0.5, as a function of the transverse 
magnetic field. The top panel corresponds to first neighbors 
(L = 1). At the factorizing field B s = ^/J y J x w 0.71.4, E 
vanishes whereas D approaches the same finite limit (|21[) V 
L, with 8 determined by (|29|) . Bottom: Discord for all sepa- 
rations L. A finite saturation limit is approached for large L 
if \B\ < B c . 

of an arbitrary pair will be of the form ([6|), where the el- 
ements ([7])- (JHJ) can be evaluated with the expressions of 
the Appendix. 

Exact results for the discord and entanglement of pairs 
are shown in Figs. |3H5]for x = Jy/Jx = 0.5 and two dif- 
ferent sizes n. For a first neighbor coupling and even n, 
the sign of J x can be changed by a local transformation 
s 2i.n ~ s 2i,n for A' = x i Vi so that both the ferromag- 
netic (J x > 0) and antiferromagnetic (J x < 0) cases at 
fixed x wm exhibit exactly the same entanglement and 
discord. They are also independent of the sign of B. 

It is immediately seen that pair entanglement and dis- 
cord exhibit significant differences for fields \B\ < B c = 
5(1 + x)Jx (the critical field of the thermodynamic limit 
n — > 00). While in this case entangleme nt p ractically 
vanishes at the factorizing field |2fj, |2lj, |23I425| 

B s = \J JyJx , (27) 
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where the chain possesses an exactly separable and de- 
generate parity-breaking ground state [24| (see Eq. pHjl ) , 
the discord remains non-zero, reaching in fact its max- 
imum in its vicinity. Besides, entanglement rapidly de- 
creases as the separation L between the spins increases, 
being nonzero for L > 2 only in the immediate vicinity 
of B s , where it is very small. In contrast, the discord de- 
creases only slightly with separation for \B\ < B c , reach- 
ing a saturation value for large L. Moreover, it is strictly 
independent of L at the factorizing field B s . 

In order to understand these results, we recall that at 
B = B s , the uniform parity-breaking separable state 



|6) 
cos 9 



\9i), \6i) =cxp[i0a i „|O i ) 



(28) 

B S /J X = VX, (29) 
where Si Z \0i) = — ^\0i), is an exact ground state (we have 
here assumed J x > 0; for J x < 0, 9, t — > (— l) l 9). It is 
obviously a state with spins fully aligned along an axis 
forming an angle 9 with the original — z axis. Due to 
parity symmetry, the partner state | — O) = P z | 6) is also 
an exact ground state at B s . They can therefore be exact 
eigenstates of H only when levels of opposite parity cross 
[24| . The ground state of the present chain undergoes ac- 
tually n/2 spin parity transitions as the field B increases 
from (reminiscent of the m/2 S z transitions of the XX 
limit [331), the last one precisely at B s . Outside these 
transitions, the energy gap between the lowest states of 
each parity is small for \B\ < B c , but nonetheless finite 
in a finite chain. 

In the immediate vicinity of B s , the exact ground state 
will therefore possess a definite parity and the correct 
side-limits at B s will be determined by the states 



|e)±|- e) 
v/2(i±(-e|e» 



(30) 



where (— 0|0) = cos™ 9 J24J. The reduced state of any 
two spins derived from (|30p will be given precisely by 
the same mixture (fT9| if the complementary overlap 
(— 0„_2|0„-2) = cos™ -2 9 is neglected. We can then im- 
mediately understand the qualitative difference between 
entanglement and discord for \B\ < B c . As B approaches 
B s , the discord between any two spins will approach a 
common finite limit for any separation L, given by Eq. 
([2"T]) with 9 obtained from ([29]) (D w 0.145 in the case 
of fig. 21 where 9 = 7r/4 at B s ). In contrast, the pair 
entanglement will vanish for B — > B s (for a negligible 
complementary overlap) as the state (fT8|) is separable. 

In Figs. [4][5] we have taken the exact ground state with 
its correct parity. The non-negligible discord between 
any two spins for \B\ < B c can then be understood in 
a similar way, as in this region the ground state can be 
seen, approximately, as a definite parity combination (|30|) 
of "mean field" states with broken-symmetry (|28|) . with 
cos 9 = B/J x , plus additional corrections. The reduced 
state of a spin pair will then be again given essentially 
by the mixture p(9) (Eq. ([T5])) plus smaller corrections, 
with the discord arising principally from p(9) (although 
corrections are non- negligible; see Fig. [5]). 



■ o.i -. 



- 0.1 - 



„- 0.1 



Q 0.1 - 




FIG. 5. (Color online) Same details as Fig.[4]for a chain with 
n = 10 spins. The common different side-limits of D at the 
factorizing field, given by Eq. (|24|) with coherence factor (|33|) . 
are now appreciable. 



Let us remark that the same reduced density ([T8j) 
arises from the statistical mixture po = i(|0 + )(0 + | + 
|0~)(0 _ |) if the overlap is discarded. p Q represents the 
T — >• + limit of the thermal state of the system at B s . 
The limit ([21]) of the discord at B s remains then also 
valid at sufficiently low T. 

On the other hand, for strong fields B ^> J x , the 
ground state is essentially the state with all spins fully 
aligned along the — z axis plus small pcrturbative cor- 
rections. As seen in Fig. [U the discord in this region is 
rather small and decreases rapidly with separation, since 
the previous superposition effects are no longer present. 
Moreover, the discord between first neighbors is very 
close to the entanglement of formation, as verified by 
a perturbative expansion: For \B\ 3> J x and L = 1 in 
the case ([26]) we obtain, setting 77 = (J x — J y )/(8B), 



D^r) 2 (-log 2 ri 2 +log 2 e-2) 
E^n 2 (- log r? 2 + log 2 e) . 



(31) 
(32) 



Thus, E is in this region slightly greater than D, as ver- 
ified in the top panel of Figs. HHH 

Results for a small chain of n — 10 spins are shown 
in Fig. [5J Although the behavior is similar to that for 
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n = 100, finite size effects become important and the 
overlap (— 0|0) can no longer be neglected. The effects 
on the discord and entanglement of the ground state par- 
ity transitions taking place for \B\ < B s are now visible, 
giving rise to small steps in these quantities. The final 
step takes place at B s , where D exhibits now a finite 
discontinuity due to the parity splitting arising from the 
coherence term (sec. IIIDj) . no longer negligible: the ac- 
tual reduced state of a spin pair derived from the states 
([50)1 is given by Eq. (|22J) with 



±cos r - 



(33) 



where the — (+) sign corresponds to the left (right) side 
at -B s , i.e., negative (positive) spin parity. The side- limits 
of D at B s are then given by Eq. (J22]) for the values ([2~5])- 
(|55f of 6 and e (leading to £>_ « 0.153, D + w 0.137 in the 
case of Fig. [5]). Small but non-zero common side-limits at 
B s of the entanglement between any two spins also arise 
24] . as determined by Eqs. ([10)) . (|25|) . In contrast with 
the discord, E q. (1251) satisfies of course the monogamy 
inequality [28l . |29| , reaching its maximum value 2/n in 
the XX limit J y — > J x (where 9 — > and C(e) — > but 
C(— e) -> 2/n, as |0~) approaches the VF-state (2i|). 

The present values of the discord at the factorizing 
field, determined by Eqs. (|2ip or in general (|24|) . are ac- 
tually much more general: Uniform chains or arrays with 
ferromagnetic ( J x 3 > 0) XY couplings of arbitrary range 
but common anisotropy \ = Jy I Jx e (0; 1) wm & l so 
exhibit a factorizing field, given for spin 1/2 by @, HH| 



(34) 



where the chain will possess again the same degenerate 
separable ground states (|28|) . Side-limits at B s will then 
be determined by the same definite parity states (f3"U)) . 
The same occurs in XYZ arrays if \ = (Jy ~ J l J)/(Jx ~ 
J**) is constant 0, with B s = s/xT,jM J * ~ * As a 
result, the ground state pair discord in all these systems 
will be finite and independent of pair separation or cou- 
pling range in the vicinity of B s , and given by Eqs. (|21l) or 
(|2"4")h with the values (|2l?)) - (j3"3"l) . Similar arguments will 
apply in the vicinity of more general factorizing fields 

An example is provided in Fig. [6j where results for a 
fully and uniformly connected XY array (LMG model 
SI), 



J*j = (1 - 8iiW(n ~ 1) . = > 



(35) 



with the same anisotropy J v / J x = 0.5, are depicted. The 
exact result can be here obtained by direct diagonaliza- 
tion, as H can be expressed in terms of the total spin 
operators = J2i s i^- The reduced pair density will 
obviously be independent of separation at any field. The 
pairwise entanglement is then small for large n, with a 
0(n _1 ) concurrence [Hj]. 

The discord is, however, non-negligible and practically 
n-independent for large n. It is verified in Fig. [S] that at 



— 0.1 




FIG. 6. (Color online) Top: Discord between spin pairs in the 
n = 100 fully connected array (a), together with the discord 
(|2ip of the state (ffSJ) at the mean field angle cos# = B/ J x (c). 
The result for first neighbor coupling at the same anisotropy 
and size n is also depicted (b), for large separation (L = n/2). 
The limits at the factorizing field B 3 are exactly coincident. 
Bottom: Same details for n = 10 spins. Here (c) depicts the 
result (1241) for the actual mixture ()22|) . which includes the 
coherence term. Side limits at B s are again coincident. 



B = B s , the same previous limit (|2T|) (n = 100) and side- 
limits ([2i ]) -([33] ) (n = 10) at B s are obtained. Moreover, 
the simple model (|18p (large n) or ([22|) (small n) for 
the reduced density of a spin pair accurately describes 
the discord in the whole region \B\ < B c = J x (and 
not just at Bs) if the mean field value cos# = B/J x is 
employed for 8, as seen in Fig. [Sj This indicates that 
the effects on D of correlations beyond the basic mean 
field description with parity restoration of the ground 
state given by Eq. (j30|) . become in this system very small, 
being negligible for large n. In contrast, such mean field 
model cannot accurately describe the pair discord (either 
for small or large separations) in the nearest-neighbor 
chain away from B s . 



IV. CONCLUSIONS 

We have first discussed in detail the discord of a mix- 
ture of aligned pairs in two different directions. While 
exhibiting no entanglement if coherence effects between 
both directions are negligible, these mixtures do exhibit 
a finite and non-negligible quantum discord if the direc- 
tions are not coincident nor opposite. In the presence of 
coherence terms, however, entanglement becomes finite 
and can bc larger than the discord. 

Such mixtures become of crucial importance for study- 
ing the pair discord in the exact (and hence of definite 
parity) ground state of finite XY and XYZ spin chains 
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in a transverse field. They represent the actual reduced 
state of an arbitrary pair in the vicinity of the factoriz- 
ing field B s . Previous results imply then a finite discord 
between any two spins in the vicinity of B s , irrespective 
of pair separation or coupling range. These mixtures are 
also the main term of the reduced pair state in the whole 
region \B\ < B Cl implying there a non- negligible pair dis- 
cord even for pairs not linked by the couplings, as was 
seen in the nearest-neighbor case. Such mixtures do in 
fact accurately describe the pair discord V \B\ < B c in 
the fully connected XY model. 

The behavior of the discord, which is free from the 
monogamy restriction, differs then considerably from 
that of the pairwise entanglement, whose limits at B s 
are small and determined by the coherence term. This 
term gives rise to a parity splitting and hence to a finite 
discontinuity of the discord at B s , visible in small chains. 

A final remark is that the present results, together with 
those previously obtained for the entanglement [HI , allow 
to identify the factorizing field as a quantum critical point 
for the ground state of an XY or XYZ chain of small size: 
At B s , the last ground state parity transition takes place 
and in its immediate vicinity, pair quantum correlations, 
as measured by the discord, become independent of both 
pair separation and coupling range. 

The authors acknowledge support of CONICET (LC, 
NC) and CIC (RR) of Argentina. 



Appendix A: Exact solution of the finite cyclic XY 
chain 



The Jordan- Wigner transformation (36[ allows to 
rewrite Eq. ([3]) for the case of first neighbor couplings 
([2"o| and for each value (±) of the spin parity P z (Eq. 
dU)), as a quadratic form in fermion operators cj, Cj de- 



fined by cj = sf exp[— wr^} = i s^s 

n 

H ± =^ B{c\ Ci - |) - |^(J +C t Ci+1 + J-c\c\ +l + h.c.) 

i=l 

= 5>*(°fc°*-5). (Al) 
keK± 

where J± = ^(J x ± J y ) and, in the cyclic case n + 1 = 1, 
Vi = 1, vt = 1 - 25 ™ M- In dSD, 

A 2 . = (B — J + cos w fe ) 2 + J 2 sin 2 u> k , uj k = 2irk/n , 

with K+ = {|, ...,n — \}, K- = {0, ...,n - 1}, 
i.e., k half-integer (integer) for positive (negative) par- 
ity [13, [37|. The diagonal form (TATT) is obtained 
through a discrete parity- dependent Fourier transform 



i— l + — 1 



k/4 



E 



feeif± 



c'u, followed by a BCS trans- 



formation c 



't 



u k a' k + v k a n - k , c n _ k 



u k a„ 



v k a\ 



to quasiparticle fermion operators a k , a k , with n k ,v k = 
|[1 ± (B - J + cosuj k )/X k }. For B > we set A fc > 
for k 7^ and Ao = J+ — B, such that the quasiparticle 
vacuum in H~ is odd and the lowest energies for each 
parity are E ± = -\Y.kek ± At B = B » = \JJxJ y , 



Afc = ,7+ — B s coswfc and E = —nJ+/2 [24 1. 

The concurrences in the fixed parity GS can be ob 
tained from the contractions /; = (c\cj)± — gi = 

J J 



{c}cj)± and the use of Wick's theorem [36|, leading 
to (4) = f , (si4) = /o 2 - ff + and at 



I l T ill alL ^ "J 

\[det(A+) T det(A-)], with (Af )y = 2(/ i _ i ±i+ 5i _ i ±i) 
I x I matrices. All previous results have been explicitly 
checked for small n with those obtained from a direct 
diagonalization. 
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